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We investigate classical scattering off a harmonically oscillating target in two spatial dimensions. 
The shape of the scatterer is assumed to have a boundary which is locally convex at any point 
and does not support the presence of any periodic orbits in the corresponding dynamics. As a 
simple example we consider the scattering of a beam of non-interacting particles off a circular hard 
scatterer. The performed analysis is focused on experimentally accessible quantities, characterizing 
the system, like the differential cross sections in the outgoing angle and velocity. Despite the absence 
of periodic orbits and their manifolds in the dynamics, we show that the cross sections acquire rich 
and multiple structure when the velocity of the particles in the beam becomes of the same order of 
magnitude as the maximum velocity of the oscillating target. The underlying dynamical pattern is 
uniquely determined by the phase of the first collision between the beam particles and the scatterer 
and possesses a universal profile, dictated by the manifolds of the parabolic orbits, which can be 
understood both qualitatively as well as quantitatively in terms of scattering off a hard wall. We 
discuss also the inverse problem concerning the possibility to extract properties of the oscillating 
^ ' target from the differential cross sections. 
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I. INTRODUCTION 



In a scattering process the interaction of the incoming projectiles with the target is of spatially local character, 
finally yielding free outgoing particles. The main question is to explore the imprints of the interaction potential on the 
K% ' outcoming beam. The basic tool for this study are scattering functions and cross sections relating physical quantities 
characterizing the outgoing beam with associated quantities of the incoming beam. The scattering angle or the escape 
(dwell) time as a function of the appropriately defined impact parameter are typical examples of scattering functions 
while the distribution of the deflection angle of the outgoing particles typically defines the angular cross section. 
Scattering problems occur in many areas of physics and chemistr y m , Hf including celestial mechanics [H, 3, [a, [y, 3, [Hi , 
charged particle trajectories in electric and magnetic fields [3, [l(|, hydrodynamical processes [Tu \1% [lj], models 
of chemical reactions [lj], 

MM 

Il7j and scattering in atomic, molecular and nuclear physics jl8l . Il9j | . When the 
scattering process evolves in two or more spatial dimensions and the interaction potential is nonlinear, the scattering 
functions may possess complex structures formed by a fractal set of singularities, originating from homoclinic and 
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heteroclinic intersections of the stable manifolds of unstable periodic orbits (UPOs) occurring in the dynamics of the 
system. The appearance of such a set of singularities is described by the term chaotic scattering [![, an active field 
of research du ring the last two decades. In general this set of singularities possesses a non-zero fractal dimension. As 
shown in Ref. [2(| the conclusion on the existence of chaotic scattering via singularities of the scattering functions can 
be a subtle issue. This is because there exist dynamical systems [2(| possessing a countable, infinite, self-similar set 
of singularities which do not exhibit topological chaos. On the level of cross sections, which are the most appropriate 
quantities to be studied in a scattering experiment, the fractal sin gula rities of the scattering functions induce a set of 
rainbow singularities (RSs) which also possess a fractal geometry [2l|, |22| . 

Most of the work done within the framework of chaotic scattering (or scattering in higher dimensions) deals with 
a static scattering potential. Recently, it has been found that complex behavior with different characteristics can 
be observed in scattering systems involving a time-dependent scattering potential [H, [24], [2^, [2(| [53, [H, [29[ . As 
a representative example of such a system the scattering of freely moving, non-interacting particles off two hard, 
infinitely heavy, oscillating discs on the plane has been studied. The absence of an energy shell in these systems 
leads to an increase of phase space dimensionality. The explicit time-dependence of this system leads to a reduction 
of the impact of the UPOs and their manifolds on the properties of the resulting scattering functions. In Ref. [25| 
the term dilute chaos was introduced to describe the behavior characterized by the accumulation of peaks, associated 
with processes leading to maximum energy exchange between scattered particles and an oscillating target around the 
position of the UPOs of the system. The set of the observed peaks possesses a self-similar structure either in an 
approximative sense occurring only between two scales determined by the geometry of the particular setup 24j or in 
a local sense in the neighborhood of specific isolated points in phase space [25| ■ The significant role of processes with 
maximum energy exchange in the scatterin g of f a time-dependent hard potential becomes more transparent in the 
case of scattering off a single oscillating disc [26| . Despite the absence of UPOs in this system the scattering functions 
possess a nontrivial structure dictated by elementary processes involving single or multiple collisions between the 
projectile and the oscillating target. Associated with these collision events are processes for which the scattered 
particle escapes from the interaction region with minimal outgoing velocity. These processes lead to peaks in the 
scattering function describing the dwell time as a function of the initial velocity of the incoming particle. The term 
low velocity peaks (LVPs) was used in Ref. [24| for the description of these structures. In phase space terminology 
the LVPs are related to the approaching of the trajectories of the scattered particles to the parabolic orbits and their 
manifolds. The parabolic orbits are the dense set of phase space points lying in the configuration space outside of 
the interaction region and having velocity equal to zero. The corresponding manifolds consist of phase space points 
which approach asymptotically the parabolic orbits. A scattering orbit approaches only marginally these manifolds 
when, due to the energy loss process, it escapes the interaction region with a very low velocity (LV). 

So far the analysis of the scattering dynamics in time-dependent scattering potentials has been restricted to the 
study of scattering functions which, although being a useful methodological tool, are not easily accessible with respect 
to experimental observation. The aim of this work is to present, for the first time, calculations of differential cross 
sections, which are observable quantities, for time-dependent scattering off hard potentials on the plane. In this case, 
in addition to the cross section associated with the distribution of the deflection angle for the outgoing particles, the 
cross section associated with the distribution of the final velocity of the scattered projectiles is of particular interest, 
since the energy of the incoming particles is not conserved during the scattering process. It will be shown that when 
the velocity of the incoming particle is of the same order of magnitude as the maximum velocity of the oscillating 
scatterer, the corresponding cross sections possess a rich structure consisting of RSs, which can be understood in detail 
both at a qualitative as well as a quantitative level using the appropriate scattering functions. Our treatment reveals 
that the key mechanism leading to these singularities is related to the interplay between low-velocity and multiple 
collision processes. In addition, it will be argued that for a single scatterer of arbitrarily shaped convex boundary the 
relevant quantity associated with the structures in the cross sections is the phase of the oscillating target at the first 
impact, which depends on the shape of the scatterer. Within this simple reasoning it is possible also to investigate the 
inverse scattering problem [3(| concerning the extraction of the characteristics of the harmonic movement of a circular 
scatterer (frequency, amplitude, direction of oscillation) as well as its size, from the profile of the cross sections and 
scattering functions. 

The paper is organized as follows: in Sec. II we present the equations which determine the dynamics of the scattering 
problem. We also discuss the role of the various parameters as well as their relevance to the subsequent analysis. 
In Sec. Ill we define the scattering functions and cross sections. In Sec. IV we use a simple Gedankenexperiment 
consisting of a suitably chosen initial beam moving parallel to the oscillation axis of a driven wall, in order to reveal 
the significance of the phase of the first collision with the moving target for the description and classification of the 
dynamics in scattering off time-dependent potentials. We introduce a three-dimensional (3D) plot relating the total 
number of collisions during the scattering process with the initial velocity of the incoming particle and the phase of 
the oscillating target at the instant of the first collision, as the basic methodological tool allowing the qualitative and 
quantitative understanding of the scattering dynamics. In Sec. V we make a comparative study of the time-dependent 
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FIG. 1: An arbitrary scatterer with convex surface is drawn in two positions: when its center is at the coordinate origin (dashed 
line) and at a displaced position where an impact takes place (solid line). Indicated is also a beam of particles as well as the 
important vectors that specify the motion of the scatterer and a projectile before and after the collision event. 

scattering off an oscillating inclined wall with the scattering off an oscillating disk in order to reveal the influence of the 
geometrical characteristics of the target on the experimentally accessible observables. In Sec. VI we discuss briefly 
the inverse scattering problem using cross sections and scattering functions to extract dynamical and geometrical 
characteristics of the scatterer. Finally, we provide in Sec. VII a summary and concluding remarks. 

II. SCATTERING SETUP 
A. Scattering dynamics 

We investigate a dynamical system which, in the general case, consists of a harmonically oscillating scatterer of 
arbitrary shape with convex boundary in two dimensions, and a beam of particles incident to the interaction region, 
as shown in Fig. [TJ The interaction region consists of all the points of the two-dimensional (2D) configuration space, 
where collisions between the particles and the target are possible. The boundary of the scatterer is impenetrable and 
its mass is assumed to be much larger than the mass of the particles. The center of the scatterer, which coincides 
with the coordinate origin when the scatterer is at its equilibrium position, moves according to the law: 

r s = Asin(wi + <f>o) (1) 

where lu is the frequency, A is the amplitude vector and 4>q is the initial phase of the oscillation. 
The oscillation angle formed between the positive x semi-axis and the scatterer's oscillation axis is: 

a = tan^(^j (2) 

The scattering dynamics of a particle with position r and velocity u is determined by the equations: 

r„+i = r„ + u n At n (3) 

u n+ i = u„ 2[n ■ (u„ - u s (i n+ i))]n, (4) 

where n — 0, 1, 2, £ — 1 is the index for the n-th collision, At n = t n+ \ — t n is the time of free flight in between 
collisions, u s (t) ~ Aw cos(uit + </>o) is the velocity of the scatterer and n is the normal to the target's boundary at the 
point of impact. Finally, £„ = ujt n + </>o denotes the phase of the scatterer upon the n-th collision. Note that n = 
corresponds to the emission time, and thus Ato = t\ — to is the time between emission and first collision. The last 
collision I before the particle escapes from the scattering region, determines the final outgoing velocity u OM t = ug. 
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The process where only one collision takes place before the particle leaves the scattering region (£= 1), is the most 
common one, but there are also processes where more than one collision occurs. In fact we can classify the dynamics 
of the system into classes of collision events (£CEs) determined by the total number I of collisions before leaving the 
scattering region. 

The equation that determines the time t n+ \ of the (n + l)-st collision is 

r„ + u„ At„ = A sm(u)t n+1 + </> ) + r bd , (5) 

where = (xbd,ybd) is the position vector from the center of the scatterer to the point of impact (see Fig. []}. 
Equation ([5]) is solved numerically, keeping the smallest positive solution. 

B. Beam configuration 

The following assumptions specify a typical scattering experiment (see also Sec. V). The beam of identical projectiles 
is synchronized with the scatterer's oscillation, i.e. the emission time to — is common for all particles and we 
choose <j) = without loss of generality. The initial position vector is r = (x ,b), where x is the common 
x-coordinate of all particles and b is the impact parameter, which defines the corresponding y-coordinate. We use 
xo <C Xbd,min — \A X \, where Xbd,min is the minimal x-coordinate of the boundary. We also assume a uniform distribution 
of&e [ 

ybd,min ~ \A y \, ybd,max + 1 A y |] , such that the entire projection of the interaction region onto the y-axis is covered, 
as ybd,max and ybd,min are the maximal and minimal y-coordinates of the boundary. All particles of the beam initially 
move parallel to the axis of oscillation, which is chosen to be the rc-axis (a = 0), and thus Uo = (uq, 0), A = (A x , 0), 
and b G [ybd,mim ybd,max]- The initial magnitude of the particle velocity uo is the same for all particles of the beam. 
This choice is compatible with the profile of monochromatic beams often used in scattering experiments since they 
facilitate the detection of the characteristics of the target. 

C. Vectors and parameters of importance 

Three vectors determine the setup of the scatterer and the beam: the initial velocity of the beam particles u , the 
amplitude vector of the oscillation A, and the normal to the target's boundary at the point of impact n. The first 
two, have been already chosen parallel to each other lying on the x-axis. We will consider in Sec. VI the case of an 
arbitrary oscillation angle a^O, in the context of the inverse problem. The normal n is related to the shape of the 
target. We intend to address the following cases: 

1. n = constant on the flat surface of an inclined wall (Sec. V A), 

2. n = n(6) on the curvilinear boundary of a disk (Sec. V B). 

We should bear in mind that in the case of a circular scatterer oscillating parallel to the beam (a = 0), n is related 
to the impact parameter b, because each particle of the beam with different b collides with the scatterer at a point 
with different n, which justifies choosing b as an independent variable to define the scattering functions. 

What is the main difference between the scattering off static and oscillating targets? When the scatterer is static, 
the outgoing velocity has always the same magnitude as the ingoing. Only the direction of the outgoing velocity is 
affected by the shape of the target, and particularly by n = n(b). When the scatterer oscillates, both the magnitude 
and the direction of the outgoing velocity depend on the target's velocity at the impact time (which in turn is 
determined by the phase of the oscillation upon collision) as well as on the shape of the scatterer and the impact 
point on its boundary (fi = n(&)). Therefore, to examine the effect of the time-dependence of a hard potential, we 
must consider orbits, which differ with respect to the first collision phase £i. We focus on £i and not on the phase of 
the second or third collision since, as will be shown in the following, £i determines to a large extent the evolution of 
the orbit and the appearance of ^CEs with I > 1. 

Following the above reasoning, it is worth determining the parameters which affect £i in the general case of a parallel 
beam of particles incident to a hard oscillating target of arbitrary shape (given w and A). These are: (a) the initial 
velocity, uq, (b) the initial phase of the scatterer when the particle is emitted, £o = Luto + 4>o, (c) the initial position 
along the ir-axis, xq and (d) the initial position along the y-axis, i.e. the impact parameter b. If one of these quantities 
changes, keeping the others constant, £i changes too. Therefore, there are four independent ways of affecting £i in 
time-dependent scattering processes. Our incoming beam is characterized by a single value uq of the velocities of all 
particles and varying uq provides us with the possibility to tune the complexity of the scattering process: Certainly the 
most interesting behavior has to be expected for velocities uq ss Auj. The remaining three parameters are equivalent 
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to each other. Let us focus on 6, which is the usually varied parameter in a typical experimental setup (see Sec. II B). 
For an oscillating scatterer of arbitrary shape, the difference with respect to £1 for each particle of the beam originates 
from the different distance which it has to travel in order to collide for the first time with the scatterer. Therefore, 
the impact parameter 6 parameterizes both n = n(6) and £1 = £1(6), and thus completely specifies the first collision 
and its subsequent dynamics. 

III. DIFFERENTIAL CROSS SECTIONS AND SCATTERING FUNCTIONS 

For the collision experiment we focus on the distributions of (a) the magnitude of the outgoing velocity u out = 

\ / u out,x u out,y 

and (b) the scattering angle 9 = tan 1 {jr^ft)- These are related to the corresponding differential 

cross sections: 



da 1 dN 

T q =<7t -N- t -^ M = Uo ^ d) (6) 



where N t is the total number of particles and at is the total cross section. Equivalently, one can calculate the 

J_ dN_ 

N t ' dq 



respective probability density functions (PDFs): g q = jj- ■ (q = u out ,6), since all the information is included in 



these quantities while at is a multiplicative factor. 

To explain the properties of the PDFs we will analyze the scattering functions: the outgoing velocity u ou t(b) and 
the scattering angle 9(b) both in terms of the impact parameter b. The following relation motivates their usefulness 
with respect to the analysis of the PDFs: 

1 dN db n . 
e « = ^-^-da- {q = U ^ 9) (?) 

where Ntr^ = \ Vbd , m J +Vbd , max is constant. 

Of particular importance for our analysis, is the appearance of smooth local extrema of the scattering functions 
at some value b = b* , such that | 6=6*= and J 1 ^"' |fc=fc*7^ 0. In this case, the corresponding PDF exhibits a 

typical square root singularity at u ou t(b*) as it can be shown by a Taylor expansion of u out (b) around b = b* up to 
order (b — b*) 2 

u ou t(b) = u out {b*) + ^± l^. -(b - b*) 2 . (8) 

If we set u ou t = u out (b), u* out — u out (b*), Au ou t — u ou t — u* out and Ab = 6 — 6*, then it follows from Eq. ([8]) that 

A b 1 

a « ; ( 9 ) 

Au out Uout - U ou t 



According to Eqs. I©-© we conclude that in the neighborhood of u* 



out 



Qu out oc (10) 

V U ut — U ou t 

Similarly, smooth extrema of 9(b) lead to square root singularities in g(9). This type of singularities in the cross 
sections (or equivalently in the PDFs), called rainbow singularities (RSs), have already been observed and analyzed in 
static scattering processes |22j | . In this case a smooth maximum of the scattering function occurs between singularities 
with fractal structure associated with the intersection of the manifolds of the UPOs of the system. The time-dependent 
setups considered here do not possess any periodic orbits. Yet, as will be shown in the following, the cross sections 
possess a set of RSs which can be attributed to parabolic orbits and their manifolds, indicating that the origin of 
these singularities is of fundamental character (see also the discussion in the introduction). Let us underline here that 
the positions 6* of the extrema of u ou t(b) and 9(b) are identical. This is an important feature that allows to connect 
certain outgoing velocities with certain scattering angles corresponding to the same RSs as it will be illustrated in 
the case of the oscillating disk (cf. Sec. V B). 
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IV. A GEDANKENEXPERIMENT TO UNDERSTAND SCATTERING OFF TIME-DEPENDENT 

TARGETS 

A. Gedankenexperiment setup 

In the present section we analyze the structures occurring in the scattering off a vertical hard wall oscillating 
horizontally according to a harmonic time-law. This setup greatly simplifies the geometry of the scatterer, yet for the 
usual beam profile of synchronously emitted particles with the same initial velocity (uq) originating from the same 
location on the x-axis (xq), it leads to a single value of £1 and consequently to trivial results: Constant scattering 
functions and delta-shaped cross sections. Since our main interest are properties in the outgoing channel associated 
with the energy of the scattered particles it seems reasonable to appropriately modify the initial beam, in order to 
allow for variations of £i. To achieve this we introduce an asynchronous particle emission while keeping both uq 
and xq constant or, equivalently, vary xq while keeping uq and £o constant. We end up with values of £1 covering 
the interval [0, 2tt) uniformly. However a beam profile containing asynchronously emitted particles or particles with 
a certain distribution of the initial distance from the target is not easily at the disposal of the experimentalist and 
should therefore be considered as a Gedankenexperiment. In the experimental setups with the inclined wall and the 
disk (see Sec. V) the interval [0, 2ir) for the values of £i is covered several times corresponding to several intervals 
of b (&-zones). We use the term 6-zone for the interval of b which corresponds to an interval [0, 2w) of £i. The 
partitioning of the initial beam into 6-zones depends on the specific geometry of the scatterer. The main purpose of 
the Gedankenexperiment is to explicitely illustrate what happens within a single &-zone. 

An important property of the vertical wall model which simplifies significantly the analysis of the scattering tra- 
jectories is that the vector n is constant and always antiparallel to the vector of the incoming velocity. Therefore, 
Eq. (HJ) becomes: 

u n +i = -u n + 2u s (t n+ i) (11) 
and the particle leaves the interaction region in opposite direction after a few collisions. 

B. Key results of the Gedankenexperiment 

We employ in the following the dimcnsionless control parameter v = -j 2 -, that is, the ratio of the initial velocity 
and the maximum scatterer's velocity. The dynamics of the scattering system is simplified in the limits vq <C 1 or 
v o 3> 1 where static scatterer approximations are valid. In particular for vq <C 1 the dynamics can be approximated 
by scattering off an effectively static potential obtained by appropriately averaging in time the exact time-dependent 
potential, while for vq ^> 1 the movement of the scatterer can be totally neglected and the corresponding observable 
quantities are similar to those obtained in the case of a static scatterer. Between these two limiting cases the complexity 
of the dynamics increases and becomes most pronounced with respect to the scattering functions as well as the PDFs 
for vq — 1. We consider the case of small oscillations of the scatterer (A — 0.01). Lengths are scaled by A, velocities 
are scaled by Alo (in which case they are denoted by v instead of u) and the frequency u) is set equal to 1. Finally, £i 
is scaled by 2tt (in the figures) with the additional convention that for phases greater than 0.5 we subtract 1 in order 
to get a symmetric representation around 0. 

The influence of the parameter t> as well as the role of £i is demonstrated in Fig. a) where the total number of 
collisions I as a function of v and £i is plotted. The results of the scattering processes with respect to the outgoing 
velocity are illustrated in Fig.[^b). In Fig.[3Ja) the outgoing velocity v ou t is presented as a function of £i € [—0.5, 0, 5) 
for Vq = 1. 

In Fig. [2|a) we observe that in the region of large Do-values only scattering processes with a single collision event 
(ICEs), as in the static case, occur. As vo decreases processes where the particle escapes after 2,3,... collisions 
gradually appear. This continues until vq — ► 1 + , where processes with an infinite number of collision events can occur. 
For vq < 1 the ^CEs with I > 1 disappear one by one rather abruptly, starting from the larger t (1 — > oo), while 
domains of inaccessible £i appear (black regions). Let us next analyze the associated scattering processes in more 
detail with respect to £i and vq: 

1. If the first collision takes place when scatterer and incident particle move in opposite directions, (head-on 
collision), i.e. for £i e [tt/2,3tt/2] (in the figures [0.25, 0.5] (J[— 0.5, —0.25]) then the particle gains energy and 
thus escapes after a single collision. This holds for any value of the parameter vq. The maximum outgoing 
velocity u out , m ax corresponds to £i — 0.5 where the magnitude of scatterer velocity is maximum | —Au}\, leading 
to u ou t t max = \ — uo — 2Aoj\. Decreasing or increasing £i starting from 0.5 the particle velocity decreases 
smoothly (see Fig. ^h) and Fig. [3Ja) for vq — 1). This simple behavior is responsible for the formation of a 
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FIG. 2: 3D plots for the vertical wall of the Gedankenexperiment: (a) The total number of collisions i as a function of vo and 
£1, (b) The outgoing velocity v ou t as a function of wo and £1. Very similar results hold for a single fe-zone of the inclined wall 
or the disk. 



smooth maximum of u ou t(b), characteristic for ICE, which in turn implies the appearance of the main RS in 
the PDF (see Fig. (He)). 

2. If the first collision takes place when the scatterer and the incident particle move in the same (positive) direction 
i.e. £i S [0, 7t/2] 1J[37t/2, 2tt] -or [—0.25,0.25] after scaling- then the particle looses energy and more complex 
processes may occur: 

(a) For a large enough value vo the particle still escapes with a single collision (I = 1). The resulting particle 
velocity is negative, i.e. the initial velocity of the particle is reversed. The wall velocity is not large enough 
to change significantly the magnitude of the particle velocity so that the oscillating wall hits the particle 
once more. 

(b) When vq decreases fCEs with £ > 1 occur as can be seen in Fig. [U(a). If the first collision takes place 
near £i = the wall velocity is large enough to lead to a significant loss of energy for the particle. For 
2CEs the particle possesses a negative velocity after the first collision and subsequently the phase of the 
wall's motion will cover the interval [0, 7r] and finally, as £2 approaches 37r/2 from below the second collision 
occurs (see also Fig. 2]) . In the same figure it is also clearly demonstrated that the 2CE domain increases 
as vo decreases. Furthermore, for smaller vo, the velocity of the target near £1 — can make the particle 
move with a positive velocity after the first collision and thus more than two collisions take place cutting 
2CEs' region (for fixed i>o) into two pieces. This division into intervals belonging to different values of £ is 
also present in Fig. HJb) . Within each interval we observe that the outgoing velocity v ou t takes on values 
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^ — 1 — i — 1 — i — 1 — i — 1 — 
0.0 0.1 0.2 0.3 

b 

FIG. 3: (a) The escape velocity v ou t as a function of £i for vo = 1 and the vertical wall, (b) The scattering function v ou t versus 
the impact parameter b for a single 6-zone for the oscillating inclined wall, (c) The same as in (b) but for the oscillating disk. 

from the maximum possible for the corresponding ZCE to the low velocity (LV) limit. For vq = 1 each 
smooth maxima of tw(£i) in Fig- EHa) corresponds to a different £CF> and we observe a global reflection 
symmetry around £i = 0. Smooth maxima belonging to neighboring arcs are separated by discontinuous 
minima (LVPs of the corresponding ^CE) . 

C. ID representation and analytical calculations 

An advantage of the simple setup considered in this section is that the dynamics can be represented graphically in 
an one-dimensional plot since it evolves exclusively along the a>axis pjlj . A representative plot of this type is shown 
in Fig. [4] where the position of the wall as a function of time t is drawn as a sinusoidal line. Particle trajectories 
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FIG. 4: ID model of the scattering off an oscillating vertical wall. The sinusoidal curve displays the position of the scatterer as 
a function of time t. The straight line segments represent the movement of different particles which are ejected asynchronously. 
The LV orbits (corresponding to ei lM (dashed line) and si,m (dashed-dotted line) near £1 = 0) which define the border between 
ICEs and 2CEs are presented. They both approach tangentially the sinusoidal curve at the second point (near phase 3n/2). 
The orbit labeled el (solid line) leads to the maximum outgoing velocity for 2CEs. 



correspond to a sequence of straight line segments. The beam can be represented in ID by employing shifts of the 
initial phase £o- 

Let us now review the above discussion of the scattering process (see section IV. B) on the basis of Fig. [4] When 
Mo is large enough the inclination is also large and, after the reflection of the initial straight line segment, no further 
intersections with the sinusoidal curve (and consequently collisions with the wall) are possible i.e. we have only ICEs. 
Moreover, it is clear that the niaxinrurn negative inclination, that is Uout^max-i 

after the collision appears at £i = n. 
Similarly the minimum negative inclination after the collision will appear when £i = 0. 

The above graphical presentation of the dynamics facilitates the derivation of analytical estimates for characteristic 
quantities determining e.g. the 2CEs region. Decreasing vq second collisions become possible. If the first collision takes 
place at £i = e\ with e\ <C 1, then the reflected orbit may possess a relatively small negative inclination (velocity) 
leading to a second collision at £2 = 37r/2 — £2, < £2 <C 1. We can approximately determine the corresponding two 
points X\ and X2, as well as the corresponding velocities of the scatterer v Sj \ and v Sj 2, by expanding the trigonometric 
functions determining these quantities in powers of Si (i = 1, 2) (up to terms of order e|): 

Xl « £1, x 2 « -1+ — , u s ,i »1 — , U s ,2 ~ S 2 - (12) 

where xt, X2 are scaled by A. 

We fix vo at a typical value allowing for second collisions to take place. Then, we seek for the phases £1 such 
that at the second point (#2) the particle trajectory is tangent to the sinusoidal curve (see Fig. 2]). There are two £1 
values, denoted as £i,^ and £\,m, fulfilling this condition. These two phases determine the lower and upper limit which 
separate the region of ICE from the the region of 2CEs for a given vo value. Additionally the particles which trace 
the border between 1 and 2CEs leave the scattering region with a low velocity (LV) close to the minimum possible 
for ICE. The following equations determine the conditions that should be fulfilled for these orbits to be realized: 

X2 — x\ = viAtiLo (13a) 
v s ,2 =vx = -v Q + 2v sA (13b) 

where Ati = (37r/2 — £2 — £i)/u>. 

Substituting x%, X2, v St \, v s _2 in Eqs. (fT3]) from Eqs. (p~2|) . we obtain a 4-th order polynomial equation for £1: 



a£i 4 + bei 3 + c£i 2 + tfei + e = 



(14) 
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where a = 1/2, b = —1, c = — 3tt/2 + vq + 2, d = vq — 1 and e = 3tt/2 + 3 + Vq/2 — (3tt/2 + 2)vq which can be solved 
analytically. We substitute for example Vq = 2.14, which allows 2CEs, and find two real solutions ex u — —0.1769 and 
£i,m = 0.4717, very close to the respective values from the numerical solution of the system being —0.1759 and 0.4709 
respectively (we use here the representation of the phases symmetrically around as done in Fig.[^a)). Having found 
the two values of ex, we determine the remaining quantities characterizing the two orbits under study such as the 
velocities ux, which correspond to ICEs with minimum outgoing velocity (LVs). Within the interval lying between 
two successive LVPs a single smooth extremum appears in the scattering function. This property of the LVPs, i.e. 
their bracketing of smooth extrema of the scattering function, makes them indispensable for understanding scattering 
processes off time-dependent hard potentials. 

It should be emphasized that immediately above the smaller solution Ex.p and just below the greater one ex t M> the 
first 2CEs occur. These orbits correspond to the LVs of the 2CEs. Note that the LVs of the 2CEs are somewhat 
larger than the LVs of ICEs and thus a discontinuity appears (see Fig. G3a)) as the almost tangential intersection 
of the curve at the second impact point will lead to a more negative inclination of the reflected final line segment, 
representing the outgoing orbit of the particle. 

We next analyze the dynamics in the region ex G ( £ i,ad £i,m) which corresponds to 2CEs. Starting just above ei !jU 
and approaching ex,M, u out initially increases, reaches a maximum (which corresponds to e± = e\ (solid line in Fig. [5]), 
£i )( u < £* < £ i,m) and then decreases to the minimum (see Fig. [3^a)) . To determine e\ we look for the maximum of 

Vout =v 2 = -vx + 2v Si2 = -(-v Q + 2v s s) + 2u s , 2 (15) 

given that 

X2 — xx = UiAtiw, (16) 
where Ati = (37r/2 — e 2 — £i)/ u ->- In other words, we look for e*, such that: 

dv n ut , 



dex 



|ei=ej " 



0. (17) 



We substitute xx, x 2 , Vs,i, v s,2 from Eqs. (TT5|) into Eqs. (|15M16p to obtain a system of 2 equations with 3 variables 
ex, £2, and v ou t- Subsequently we find the relation v out — f(ex) after the substitution of e-i. Applying the condition 
(jTTJ) yields e\ — 0.1472 for vq = 2.14, which is close to the numerical value 0.1320. This smooth maximum of 2CEs, 
i.e U2 which corresponds to £1 = e*, causes the secondary RSs in the PDFs. We remind the reader of the fact that 
the main smooth maxima that appear in the scattering functions (and correspond to the main RSs) in the PDFs are 
those of the ICEs, because they always cover a greater width with respect to £1 as they always include half of the 
complete interval of £1, namely [tt/2, 3tt/2]. Additionally, if we are interested in the secondary RSs we should study 
all smooth maxima of ^CEs with £ > 1 which appear, as it will be shown, close to those of the 2CEs. 

The last calculation concerns the critical v$ r which is the threshold for 2CEs (see Fig. [D(a)). Note that vq is a 
parameter in Eq. (jT4j) for e\. 2CEs emerge first when the two real solutions ex llh and £\ m are equal (see Fig[4j. 
Therefore we ask for the value of the parameter vq («o r )> sucn that we encounter one double real solution (ef) of 
Eq. (fT4"|) . We finally find i>g r = 2.237 which lies very close to the value 2.235 obtained numerically. 

Using similar arguments we can perform the corresponding calculations for the remainder of the £CEs (I > 2), 
which appear gradually as v$ decreases. The limits of the region of 3CEs for a certain vq can be approximated if 
we suppose that the orbits we seek have £1 near 0, second collision with a phase close to 7r/2 and a third collision 
tangential to the curve at the point with phase tt. Note that the LVs of the 3CEs are slightly greater than those of 
the 2CEs, whereas their corresponding maxima are a bit lower in close analogy with the comparison of ICEs and 
2CEs. If we proceed with decreasing further Vq, we should observe the appearance of 4CEs that will divide the area 
of 3CEs into two pieces. For the approximation of the limit orbits for £CEs with £ > 3 the same two last conditions 
with 3CEs should be satisfied, but one should further demand that more collisions occur in the phase interval [0, tt/2]. 
This process continues until Do = 1 where we can have an infinite number of collisions, i.e. the particle follows the 
orbit of the scatterer. 

For vq < 1 the crucial difference is that uq is smaller than the maximum velocity Ato of the scatterer, i.e. some values 
of £1 (especially those near £1 = 0) are no more accessible (see Fig. Ufa)). This fact causes an abrupt disappearance 
of £CEs with t > 1, starting from larger £ and proceeding to the smaller ones, so that finally the area of ICEs shrinks 
until it becomes vanishingly small as vq — > 0. 

Lets conclude our Gedankenexperiment. The value of the parameter vq which corresponds to the richest behavior 
of the scattering processes is vo — 1. Several ^CEs appear, separated from the respective LVPs, with smooth maxima 
in between (see Fig. [3](a)), which in turn correspond to RSs in the PDF g Uout (see Fig.[5jc) for example). Given that 
LVs are defined as the lowest possible outgoing velocities for an £CE scattering process, they constitute the closest 
trajectories to the parabolic orbits and their manifolds. This property of the LVPs allows us to connect the RSs 
observed in the PDF with the parabolic orbits as well as their manifolds in phase space. 
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FIG. 5: Characteristic quantities for the scattering off the oscillating inclined wall (7 = 7r/3, vo = 1): The scattering functions 
(a) iw(6), (b) 9(b) and the PDFs (c) g Vout , (d) g g . 



V. THE INFLUENCE OF THE GEOMETRY ON THE SCATTERING DYNAMICS 

A. Oscillating inclined wall 

The inclined wall has finite size (this choice simulates the finite size effect in the case of an oscillating disk) such 
that the scattering region in the y-direction is entirely covered by a beam with b G [—1,1]. We examine here the 
influence of the inclination of the wall 7 i.e the acute angle between the wall and the positive horizontal semi-axis 
x. For all particles the normal is the same n = (—sin 7, cos 7) and the unitary tangential vector is also constant 
t = (sin7,cos7). In the case of a static inclined wall we have simply u out — u n and 9 = 27, i.e. the scattering 
functions arc constant and PDFs delta function like. 
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In the case of an oscillating inclined wall (7 = 7r/3) there is a continuous variation of £1 with the impact parameter 
6. This originates from the difference with respect to the distance which the particles have to travel due to the 
inclination of the wall. Therefore, it is expected that the structure found in the Gedankenexperiment (see Fig.[3ja)), 
discussed in the previous section, appears here too. Indeed, in the scattering function of the outgoing velocity v ou t(b) 
(see Fig. Eta)) we observe a repetition of the structures belonging to a single b-zone. This is due to the fact that £1 
is not restricted to the range of one cycle [0, 2ir) (as in the Gedankenexperiment), but, depending on the inclination 
and the length of the wall, the particles can collide during several cycles of the oscillation. We can divide the impact 
parameter in 6-zones each of which corresponds to values of £1 covering -continuously but not uniformely- a period of 
oscillation. The number of the repetitions of the structure is equal to the number of the 6-zones, N zn = tan ^ U0 'ijL ■ 
A typical characteristic of the inclined wall is that 6-zones have equal length and for this reason we observe a pattern 
consisting of a repetition of a single structural element. However, this is not an exact repetition, because the finite size 
of the scatterer introduces a slight difference between the boundary zones and the internal ones. The total number 
of collisions £ as a function of Vq and £1 for a 6-zone of the inclined wall (6 S [0, tan(j)uQ^j] is very similar to the 
corresponding one for the vertical wall (see Fig.^a)). Consequently, for a single 6-zone and for vq = 1, we expect a 
similar structure in the scattering function v out (b) as that of Fig. [3ta) of the Gedankenexperiment and this is indeed 
the case as we can see in Fig. [3jb) . 

The difference between the direction of the normal n and the initial velocity uo, affects scattering angles, and thus 
7T in contrast to the case of the vertical wall. The scattering function 6(b) is shown in Fig. E^b) where we can 
observe the 6-zone pattern discussed above. 

To ease the interpretation of the result of the scattering process we can rewrite Eq. ((4]) in terms of the normal and 
tangential components: 

u„+i = (u„ ■ t)t - [n ■ (u„ - 2u a (t n+ i))]n (18) 

In case we are interested only in ICE we can put n = and u\ = u ou t- Obviously, the tangential component 
does not change whereas the normal one depends on the velocity of the wall at the impact time (£1). When moving 
initially with a velocity of opposite sign compared to that of the wall i.e. £1 € (tt/2,3tt/2) the particle gains energy 
and approaches n, so it has a larger outgoing velocity and scattering angle. Therefore, 8(b) and v out (b) possess a 
similar appearance. The maximum of v ou t and 9 corresponds to £1 = 7r, where the wall velocity is minimum. 

Vout,max = V 2 + v o(% + V ) - 2(1 + V ) COS 27 (19) 



max = tan- ( ^ + 7 )S 7% ) < ,/2 + 7 (20) 
\— 1 + (1 + v ) cos 27 J 

In Figs. [5jc) and [5jd) the PDFs Q Voui and gg are presented, respectively. We can verify that the main RS of 
the PDFs corresponds to the smooth maxima for ICE derived above. For 7 = ir/3 we obtain from Eqs. (|19[) - (j2"0|) 
Vout,max = 2.65 and 9 max /ir — 0.77 which coincide with the numerical results. It should be noted that for all the 
6-zones in the scattering functions these smooth maxima are the same, so they will contribute to the same main peak 
in the cross sections increasing its height. 

Secondary RSs correspond to the maxima of £CEs with I > 1. In g Vaut we observe secondary RSs close to v out = 2.0 
(corresponding to 2CEs maximum), close to 1.35 (3CEs maximum), and an accumulation of overlapping maxima for 
Vout £ (0.7, 1.3) but still leading to very minor peaks (£CEs maximum with t > 3). The remote peak at v out = 0.5 
corresponds to the LVP for ICE. Secondary RSs for gg are very close to the dominant ICE peak (as we can observe 
in the corresponding scattering function 0(b)) and thus contribute to the increase of the width of this peak. 



B. Oscillating disk 

Let us now examine the influence of the curvature on the structure of scattering functions and consequently on 
the PDFs by investigating the scattering off an oscillating disk (radius R = 1 and center placed at the origin). The 
normal unit vector n at each point of the scatterer's boundary depends for a = only on the impact parameter 6: 



6 2 6 . . 
n=H/l-^). (2D 



In addition the unit tangent vector reads: 
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In the case of a static disk the magnitude of the outgoing velocity is equal to the initial one (v ou t(b) = uq and 
Qua^t = 3{ u out~uo)) but the scattering angle varies smoothly because of the curvilinear boundary of the disk (variation 
of n): 



I sin (0/2) | 

Qb = 1 (24) 

This form of 6(b) remains as a backbone in the case of an oscillating disk. 

When the disk oscillates, the variation of £i, arises due to the different distances that the particles have to travel 
because of the curvature. Thus, we expect a repetition of the basic structure of the 6-zone of the Gedankenexperiment 
to occur in the scattering function v ou t(b) (see Fig. EJa)). In Fig. [3fc) the scattering function v ou t(b) for a 6-zone is 
presented, exhibiting a structure similar to that observed in Fig.[3]Ja) for the Gedankenexperiment. Nevertheless, each 
particle collides at a different point on the scatterer with a different normal and tangential vector and not a constant 
one as is the case for the wall. This fact influences the values of v ou t for each 6-zone, as well as the partitioning into 
6-zones, which are not anymore of the same length, but become narrower as we approach b — R or b = —R. Namely, 
the upper limit b z of the z-th zone (z — 1, 2, .., N zn = 2 {L ) is: 



b z = \lzu Q — (2R-zuQ—) (25) 

Using Eq. ([23)) we can calculate the length of each zone: b z — b z -\ with &o = 0. For the zones b < the same 
relation (r.h.s of Eq. (|25l multiplied by —1) holds for the lower limit of the zone. Equation ([25]) also describes the 
accumulation of the maxima towards the boundaries of b (— > ±R). 

In Fig. EJb) the scattering function 6(b) is presented. One can observe smooth minima and maxima for b < and 
b > 0, respectively. This happens because n points to different directions below and above 6 = 0, namely towards 
the scattering angles (it, 37r/2) and (ir/2,n), respectively. The latter affects the scattering angle of a particle that 
(for example) approaches n (gains energy) in the following way: for b > the scattering angle becomes larger as 
the trajectory after the collision approaches n, whereas for b < it becomes smaller. The inverse happens when the 
particle diverges from n because of energy loss. 

We stress once more that there is a one-to-one correspondence of the structures (location of extrema and LVPs) 
appearing in 6(b) and v out (b). Therefore the RSs of the PDFs g Vout and gg (Fig.[6jc) and[(Jd), respectively) originating 
from the same extrema of the corresponding scattering functions are uniquely related to each other. Another important 
characteristic of the PDFs g Vmtt and ge in the case of the oscillating disk, is that the number of the main RSs 
which correspond to the smooth extremum of each 6-zone for ICEs, are equal to the number of the b-zones. This 
discrimination of RSs, which did not exist for the inclined wall, originates from the curvature (change of n) which 
forms unequal values of the extrema in the scattering functions for each 6-zone (see Fig. Uta)). Besides the extrema 
originating from ICEs, which correspond to the main RS, there occur extrema within each 6-zone in the scattering 
functions which stem from ^CEs with I > 1. These extrema have different but very close values in each 6-zone (see 
Fig. HJa)), and thus contribute to the wide peak of g Voul at approximately v ou t — 1.0 (see Fig.JBJc)). The fluctuations 
in this diagram for < v out < 0.8 stem from events with minimum outgoing velocity (LVPs). We can verify this by 
considering the corresponding scattering function v out (b) shown in Fig.[6ja). The analysis of ge (Fig.[6jd)) is similar 
to the preceding one for g Uout ■ We observe RSs corresponding to the maxima and minima of the scattering function 
with a symmetry around 6 = tt. Each peak of this diagram corresponds to a peak of g Vout and thus in the case of an 
oscillating disk, we can redistribute the particles of a plane wave like beam to certain angles and associated velocities 
within the scattering process. 

An alternative way of interpreting the scattering functions and PDFs of the oscillating disk is presented in the 
Appendix. If we are interested only in the main peaks of the PDFs resulting from ICEs, then we can use the 
approximation that the disk does not move in configuration space but only changes velocity according to the law 
Ud = Au cos(ujt + <f) ) (static disk approximation SDA) [H, HH . We derive the following analytical expressions for u out 

(l/R)(R 2 ul + 2Ajcos(C) x 
(2(b 2 - R 2 )u cos(a) + 2bVR 2 - b 2 u sin(a) + Alu cos(C) x 

(R 2 + (R 2 - 2b 2 ) cos(2a) - 2bVR 2 - b 2 sin(2a)))) 1 / 2 , (26) 
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FIG. 6: Characteristic quantities for the scattering off the oscillating disk (a = 0, vo = 1): The scattering functions (a) v ou t{b), 
(b) 0(b) and the PDFs (c) Q Vout , (d) g g . 



where C = and 



and 



6i(b) = tan,- 1 



02(b) = tan 



2b[u a VR 2 - b 2 + Amcos(()(bsin(a) - cos(a)VR 2 - b 2 )} \ 
(2b 2 - R 2 )u - 2Aujcos(0(bsm(a)VR 2 - b 2 + cos(a)(b 2 - R 2 )) J 

' 2A^cos(C)(cos(a)y / # ! - b 2 - frsin(cv)) \ 

R 2 u + 2Ajcos(C) cos(a)(b 2 - R 2 ) + 2Aujbcos(() sin(a)VR 2 - b 2 ) ' 



(27) 



(28) 



for g > and g < respectively, g = (2ud — uo) • n. Please note that cos(£) appears both in u ou t(b) and 0(b). If we 
are interested in the values of the impact parameter b — b* where an extremum of the scattering function is located, 
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we require dc ° s ^ \h=h'= that is: 



db 



b*u> sin(C) 



= 



uoV-K" - b* 2 

An obvious solution is b* = 0. The other solutions can be derived from the condition £ = rmt which gives: 



b* = ±\ R 2 



to 2 7t 2 Uq 2m7TMQ0o u o^o 2miruoXo 2uqXq(/)o 



(29) 



(30) 



where m is an integer between 



uq (ftp — Rlu-\-xqlu 

7TUQ 



and 



Ugcj> +R(Aj + XoLJ 
7TU() 

In Fig. [7] we present the values given by the exact numerical solution and the SDA for b* and v out (b*). We observe 
an overall good agreement. These values of b* and v ou t(b*) are only approximately independent of £o- However, it 
should be stressed that for values of u ou t(b) close to b = the SDA ceases to be a good approximation. This is because 
the effect on v out of the disk's motion in configuration space (which is neglected by SDA) becomes comparable with 
the effect of the curvature. For the same reason the dependence on £o becomes very important for these values. On 
the contrary, near b — *■ R the curvature forms to a large extent the outgoing velocity of the particle and SDA predicts 
accurately u out , which is also approximately independent of £o- 



VI. THE INVERSE SCATTERING PROBLEM 



Here we assume that the shape of the hard scatterer is circular and the oscillation law is harmonic. An experi- 
mentalist could estimate uj using a pulsed beam with certain repetitions in time. If the corresponding time intervals 
become equal to the period of the oscillation the synchronization will result in discriminating the peaks of the PDFs 
that otherwise should have been smoothened out due to the phase averaging. The oscillation amplitude A can be 
obtained from the maximum outgoing velocity observed, which always obeys the equation u ou t. ma x = uq + IAiq. The 
radius of the scatterer R is related to the number of dominant peaks in the PDF for v ou t > 1 (see Fig. [HIc),(d)) which 
is equal to the number of &-zones provided by the equation N zn = m 2 ^. . 

As soon as w, A, and R have been obtained there is only one information missing to complete the description of the 
disk's motion, namely the oscillation angle a. If we encounter a nonzero angle a we observe a breaking of the reflection 
symmetry around b = in the scattering function v out {b) (see Fig. |S{a) for a — 7r/3) with corresponding results for 
the PDF Q Vout (Fig. EJc)). We can quantify this asymmetry by the position b a of the point which is determined via 
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FIG. 8: Characteristic quantities for the scattering off the disk with oscillation angle a = 7r/3: The scattering functions (a) 
v out {b), (b) 6(b) and the PDFs (c) g Vout , (d) g e . 

the condition u s • n = leading to zero momentum change of the projectile due to the motion of the disk (see Eq. HJ. 
This yields: 

b a = i?cosa. (31) 

For a — 0, bo — R corresponds to the accumulation of the RSs (see Fig. [6ja) ) . For a ^ 0, b a provides us with the 
position of the minimum of the upper envelope of v ou t(b) moving closer to b = with increasing a (see Fig. 0a)). 
We stress that the variation of the normal vector n with varying b is greater in the vicinity of b = R compared to the 
neighborhood of b = 0. As a result the shift of the accumulation point b a towards 6 = affects the overall appearance 
of the scattering function which exhibits smooth oscillations near this point. This is because near the point on the 
surface corresponding to b a the effective velocity of the disk, i.e. n • u S Tl+1 in Eq. (0| is very small and the dynamics 
approach the scattering from a static disk. Consequently the region of smooth behavior of e.g. v ou t(b) near b a is 
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shifted closer to b — and its size increases. 

Regarding the PDF the variation of the oscillation axis of the scatterer, does lead to changes of both the width 
and the location of the RSs. However, the location of these peaks also depends on the initial phase of oscillation £o- 
This behavior becomes even more evident as b* approaches (see Sec. V B). Moreover a phase averaging leads to 
a smooth PDF without dominant peaks. Therefore, the location of the RSs with respect to Q Vout does not possess a 
one-to-one correspondence with the oscillation angle, and consequently cannot be used for determine it, at least for 
the case of an unknown £o- 

The effect of the rotation of the oscillation axis is shown in Fig. [5£b) for 6(b) where we observe a shift of the position 
b a of the minimum of the upper envelope, and in Fig[8jd) for the PDF where a breaking of the symmetry with respect 
to — 7r occurs. 



VII. CONCLUDING REMARKS 

In the present work, we have investigated the scattering of a beam of non-interacting particles off an oscillating 
target in the plane. Our study focuses on the exploration of the basic scattering mechanisms and on the experimentally 
accessible differential cross sections for the outgoing velocity of the scattered particles and the deflection angle. At 
the methodological level we have elevated the importance of the phase of the oscillating target at the instant of the 
first collision for the description and classification of the collision events comprising the scattering process. In order to 
gain additional insight into the scattering dynamics we have considered the simplest possible system consisting of an 
oscillating wall and a beam of particles moving with a constant velocity parallel to the axis of oscillation of the wall. 
We have shown that within this simple example the complete complexity of the scattering dynamics is recovered in a 
Gedankenexperiment with asynchronous ejection of the particles in a beam. In this system the ensemble of the particle 
trajectories is divided in subsets each containing orbits with a fixed number of total collisions between particle and 
oscillating wall. This classification has been clearly illustrated in a 3D plot displaying the total number of collisions I 
between particle and wall as a function of the initial velocity Uq of the particle and the phase of the scatterer oscillation 
at the instant of the first collision £i. A closer investigation of the underlying dynamics reveals a universal behavior 
expressed through specific critical points separating regions in the (uo,£i) plane formed by trajectories with different 
values of I. This behavior turns out to be generic for scattering off oscillating targets as suggested by our analysis 
of the scattering dynamics for two other examples, the inclined wall and the disk. In addition we have shown that 
each zone with a fixed number of collisions and £i € [0, 2ir] is characterized by a smooth maximum in the scattering 
functions relating the impact parameter of the incoming beam of particles with the outgoing velocity or the deflection 
angle. Each smooth maximum leads to a RS expressed through a prominent peak with prescribed width in the 
associated differential cross section (or equivalently the PDF). The position of the peak in the cross section is given 
by the phase space location of the maximum while its width is determined through the locations of the events leading 
to minimum outgoing velocity, bracketing the smooth maximum, in the related scattering function. It is the latter 
property which allows us to connect the peaks occurring in the cross section of the considered scattering processes 
with the parabolic orbits and their manifolds in the corresponding phase space. In fact the region of multiple collisions 
is accessible when the velocity of the beam particles becomes comparable with the maximum velocity of the scatterer. 
In this case the cross sections of the scattering process attain a characteristic profile, generated by the superposition of 
peaks attributed to RSs. The complexity of the observed pattern increases as the velocity of the particles approaches 
the maximum velocity of the oscillating scatterer due to the accumulation of the peaks in this limit. It is a unique 
characteristic of our time-dependent system that the RSs in the cross section of the outgoing velocity possess a one- 
to-one correspondence with those of the scattering angle. It is important to notice here that both the RSs in the cross 
section, as well as the singularities in the scattering functions form a finite set. No topological chaos is present in the 
scattering off the harmonically oscillating disk. 

Our analysis is valid for hard scattering off an arbitrarily shaped oscillating target with convex geometry. However, 
for a target of general convex shape the interval [0, 27r) for the phase £i can be covered several times depending also on 
the profile of the incident beam. If different profiles of the beam, with certain phase or velocity distributions are used, 
it is interesting to investigate the imprints that time-dependence leaves on the cross sections. The understanding 
gained by the present study opens up the perspective to conduct the inverse scattering problem for scattering off 
time-dependent targets. As shown this can be achieved by manipulation of the incident beam. With the appropriate 
tuning it is possible to extract information concerning the frequency, the amplitude and the axis of oscillation as well 
as the size of the scatterer. It is also interesting to study if the oscillation law, assumed to be harmonic in all the 
systems considered here, leaves its imprints in the peaks of the differential cross sections. Finally, it is of relevance to 
extend the present study to different setups supporting also the presence of UPOs in the scattering dynamics. 



18 



ACKNOWLEDGMENT 



The authors thank Dr. P. K. Papachristou for helpful discussions. L.B. acknowledges financial support from the 
projects IN-111607 (DGAPA-UNAM) and 79988 (CONACyT). P.S. acknowledges financial support by the Deutsche 
Forschungsgeeinschaft under the contract Schm 885/13. 



[1] E. Ott and T. Tel, Chaos 3, 417 (1993). 

[2] E. Ott, Chaos in dynamical systems, second ed., Cambridge University Press, Cambridge, 2002. 
[3] J. Hietarinta and S. Mikkola, Chaos 3, 183 (1993). 

[4] J.-M. Petit and M. Henon, Astronomy and Astrophysics 188, 198 (1987). 
[5] P.T. Boyd and S.L.W. McMillan, Chaos 3, 507 (1993). 

[6] J.K. Moser and P.J. Holmes, Stable and random motions in dynamical systems: With special emphasis on celestial me- 
chanics, Princeton University Press, 2001. 
[7] L. Benet, D. Trautmann and T.H. Seligman, Celest. Mech. Dyn. Astron. 66, 203 (1997). 
[8] L. Benet, T.H. Seligman and D. Trautmann, Celest. Mech. Dyn. Astron. 71, 167 (1998). 
[9] A. A. Chernikov and G. Schmidt, Chaos 3, 525 (1993). 
[10] W. Breymann, Z. Kovacs and T. Tel, Phys. Rev. E 50, 1994 (1994). 

[11] E.A. Novikov and LB. Sedov, Zhurnal Eksperimental noi i Teoreticheskoi Fiziki 77, 588 (1979). 
[12] H. Aref, J. Fluid Mech. 192, 115 (1984). 

[13] V. Rom-Kedar, A. Leonard and S. Wiggins, J. Fluid Mech. 214, 347 (1990). 

[14] C.C. Rankin and W.H. Miller, J. Chem. Phys. 55, 3150 (1971). 

[15] D.W. Noid, S.K. Gray and S.A. Rice, J. Chem. Phys. 84, 2649 (1986). 

[16] P. Gaspard and S.A. Rice, J. Chem. Phys. 90, 2225 (1989). 

[17] B.-P. Koch and B. Bruhn, Chaos 3, 443 (1993). 

[18] D. Wintgen, K. Richter and G. Tanner, Chaos 2, 19 (1992). 

[19] J.-M. Yuan and Y. Gu, Chaos 3, 569 (1993). 

[20] C. Jung, C. Mejia-Monasterio and T. H. Seligman, Phys. Lett. A 198, 306 (1995) 
[21] C. Jung and H.J. Scholz, J. Phys. A: Math. Gen. 20, 3607 (1987). 

[22] C. Jung, G. Orrelana Rivadeneyra and G. A. Luna-Acosta, J. Phys. A: Math. Gen. 38, 567 (2005). 
[23] A. Antillon, J.V. Jose, and T.H. Seligman, Phys. Rev. E 58, 1780 (1998). 

[24] P.K. Papachristou, F.K. Diakonos, E. Mavrommatis and V. Constantoudis, Phys. Rev. E 64, 016205 (2001). 

[25] P.K. Papachristou, F.K. Diakonos, V. Constantoudis, P. Schmelcher and L. Benet, Phys. Rev. E 70, 056215 (2004). 

[26] P.K. Papachristou, F.K. Diakonos, V. Constantoudis, P. Schmelcher and L. Benet, Phys. Lett. A 306, 116 (2002). 

[27] L. Benet, J. Broch, O. Merlo and T.H. Seligman, Phys. Rev. E 71 036225 (2005). 

[28] O. Merlo and L. Benet, Celest. Mech. Dyn. Astron 97, 49 (2007). 

[29] L. Benet and O. Merlo, Phys. Rev. Lett. 100, 014102 (2008). 

[30] C. Jung, C. Lipp and T. H. Seligman, Ann. Phys. (NY.) 275, 151 (1999). 

[31] Y. C. Lai and C. Grebogi, Int. J. Bifurcation Chaos Appl. Sci. Eng. 1, 667 (1991) 

[32] A. J. Lichtenberg, M. A. Lieberman, Regular and Chaotic Dynamics, Applied Mathematical Sciences 38, Springer Verlag, 
New York, 1992. 

[33] M. A. Lieberman and A. J. Lichtenberg, Phys. Rev. A 5, 1852 (1972). 



